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Abstract 

Using the continued-fraction method we solve the Caldeira-Leggett master equation in the phase-space ( Wigner) represen- 
tation to study Quantum ratchets. Broken spatial symmetry, irreversibility and periodic forcing allows for a net current in 
these systems. We calculate this current as a function of the force under adiabatic conditions. Starting from the classical 
limit we make the system quantal. In the quantum regime tunnel events and over-barrier wave reflection phenomena modify 
the classical result. Finally, using the phase-space formalism we give some insights about the decoherence in these systems. 
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1. INTRODUCTION 

Transport properties in periodical structures with 
broken spatial symmetry (ratchet systems) have at- 
tracted a lot of attention in the past decade [1] . These 
systems, under out of equilibrium conditions, can dis- 
play net current even when the time and space average 
of all applied forces is zero. For this reason, these sys- 
tems are also called Brownian motors [2] . Realizations 
of Brownian motors can be found, for example, in bio- 
logical and condensed matter systems [3] . 

In the classical limit, the dynamics is governed by 
classical Langevm or equivalently Fokker-Planck equa- 
tions. This limit has been extensively studied in the 
recent past, establishing the main phenomenologies 
of the ratchet effect (non vanishing rectified current, 
etc..) [1]. The inclusion of quantum fluctuations in 
the system enriches the ratchet phenomenology [4]. 
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However, the difficulties to deal with quantum dissi- 
pative systems, make their study harder and one finds 
less works in this regime [3,4,5,6,7,8,9]. 

Recently, it has been shown that it is possible to 
solve Quantun Master Equations using continued 
fraction methods [10,11]. Such Master Equations are 
the "quantum generalization" of the classical Fokker- 
Planck equations for some ranges of the parameters 
(weak coupling, high temperature,...). Adapting this 
technique from the classical case [10], the authors of 
[11] have studied the ratchet effect as a function of the 
temperature. In this work we extend the study to the 
force dependence of the rectified velocities in the high 
temperature/weak coupling regime. 



2. THEORETICAL FRAMEWORK 

The out of equilibrium dynamics of classical sys- 
tems is well established since one century ago with the 
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seminal works of Einstein and Langevin. The Langevin 
(Fokker- Planck) equations describe satisfactorily the 
stochastic behavior of these systems. However, their 
quantization represents in most cases a difficult task 
[12]. The most satisfactory approach consists of quan- 
tizing the Hamiltonian of the System plus its Environ- 
ment. The minimal Bath model is a large collection 
of harmonic oscillators. This approach, in the classical 
limit, recovers the Langevin description of open sys- 
tems. 

In the quantum realm a closed evolution equation 
is not possible in the general case. This is mainly due 
the non- Markovian nature of the quantum correlations. 
Nevertheless, under some approximations, it is possi- 
ble to derive a Markovian Quantum Master Equation. 
Thus, it is typically assumed that the relaxation times 
are much longer that the quantum correlations of the 
bath, 

1 h 

7 ^ T' « 
with 7 the damping parameter measuring the strength 
of the coupling to the bath. Under these conditions 
(high temperatures and/or weak coupling) the evolu- 
tion is given by the celebrated Caldeira-Leggett master 
equation [13], 



ihdtQ = [H s , g] + -^[x, {p, £>}] 



i'ymk^T 



c,g]]. 



(2) 

Here q is the reduced density matrix, obtained af- 
ter tracing out the environmental degrees of freedom 
and Hs is the system Hamiltonian, in our case Hs = 
p 2 /2m + V(x). The first term in (2) gives the unitary 
evolution, the second yields the dissipation and the 
last is the responsible for the diffusion. 

An alternative description of quantum systems is 
provided by the phase-space ( Wigner) formalism. The 
central object is the Wigner function, denned in terms 
of the density matrix as [14] 

W ^P) = i jdye- ip y^g(x + y/2,x~y/2). (3) 

This representation not only gives a phase-space de- 
scription for quantum systems, but it allows to extend 
classical concepts and tools into the quantum domain. 
Besides, it permits to find quantum analogues to classi- 
cal phenomena (Quantum Chaos, phase-space trajec- 
tories, ...) and provides a natural Quantum-Classical 
connection. 



In the Wigner representation the master equation 
(2), can be written as 

dtW =[- P d x + V d P + 7 t d v (p + dp) 

+ ^=0 ](i ( S ) y (2 s+ l) wa (2 !+ l)j ? (4) 

The first two terms correspond to the Poisson bracket 
(Liouville equation). This Poisson bracket plus the 
third term yields the classical Fokker-Planck equa- 
tion, evidencing the classical character of the noise 
and relaxation. The last term is purely quantum and 
comes from the unitary evolution of the closed system 
( Wigner-Moyal bracket). 

Equation (4) is written in dimensionless form [11]. 
The variables have been scaled with help from a refer- 
ence length xo, mass M and frequency fio- For exam- 
ple, action variables are scaled by the characteristic ac- 
tion So = MfloXQ, energy by E = Mfigx 2 , etc. Then 
the coefficient in the quantum sum is given in terms of 
the de Broglie wave length by, 

M 



K W = (-l)'A&/(2s + l)! 



(5) 



Besides, ft is introduced in terms of So via the quantum 
parameter K: 



ft/So = 2ir/K and A dB cx l/K 



(6) 



Note that the classical limit is naturally recovered let- 
ting K — > co. 

To conclude, the calculation of observables in the 
Wigner formalism is an example of extension of classi- 
cal methods to the quantum. The Quantum Mechani- 
cal expectation value is obtained via W(x,p) as a clas- 
sical "average" 



(A) =Tt(qA) = JdxdpW(x,p)A( 



x,p) , (7) 



with A(x,p) being the classical observable correspond- 
ing to the operator A (via Weyl's rule). 



3. CONTINUED FRACTION APPROACH 

A suitable non-perturbative technique to solve clas- 
sical Fokker-Planck equations of systems with few 
variables is the continued-fraction method [15]. This 
is a special case of the expansion into complete sets 
(Grad's) method to solve kinetic equations in sta- 
tistical mechanics. The technique had already been 
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adapted to quantum dissipative systems in studies 
of spins and quantum nonlinear optics. Recently this 
method has been extended to quantum Brownian 
Motion problems described by Caldeira-Leggett type 
equations exploiting the Fokker-Planck-like structure 
of the quantum master equation in the Wigner repre- 
sentation [10,11]. 

The idea of the method consists of expanding the 
desired solution W of the differential equation (d t W = 
CW) into an appropriate basis. The equations for the 
expansion coefficients (d) have then the form of a sys- 
tem of coupled differential equations, say 



Ci 



5,3+i Cj + 



(8) 



The goal is to find a basis in which the range of index 
coupling I of the coefficients C; is as short as possible. 
Indeed, for finite coupling range (/ < oo) the differen- 
tial recurrence relation (8) can be solved by iterating 
a simple algorithm, the structure of which is like that 
of a continued fraction, 



C = 



Pi 



qi + 



P2 



(9) 



92 + 
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In the case of the Caldeira-Leggett equation the ex- 
pansion is on an x & p basis: 



W(x,p) oc ^ C"u a {x)ip n {p) 



(10) 



Then, the recurrence (8) is replaced by a two index one, 
C" = 53 m p Qa™Crnj which can be transformed into a 
one index recurrence introducing appropriate vectors 
and matrices: 



(ii) 



This recurrence relation can be solved with matrix 
continued- fraction methods. In [11] the explicit form 
of matrices Q was derived. In particular, for periodic 
potentials [^(a?) = V(x + L)] explicit recurrence rela- 
tions were constructed using Hermite functions for the 
momentum basis and plane waves (u a oc e %ax ) for the 
position basis. This choice provides a finite coupling 
range /, making possible its solution by the aboves- 
ketched method. 



Solving the Quantum Master Equation in this way, 
the Wigner function is obtained [i.e., the density ma- 
trix), so any observable can be calculated. 



4. RESULTS FOR QUANTUM RATCHETS 

In the context of Brownian motion the simplest 
model that allows for broken spatial symmetry is the 
two harmonic potential, 

V{x) = -Vo[sini + (r/2) sin(2a:)] , (12) 

(another ratchet systems can be modelised using tight 
binding models, rotational spin systems, optical lat- 
tices, etc.). Potential (12) has a steeper side to the 
left (see fig. 2) , which we will call the "hard" direction, 
while we will refer to the other direction as the "easy" 
one (to the right). 

Using the continued-fraction technique we can solve 
the Caldeira-Leggett equation. In particular for the case 
of the two harmonic potential (12) the coupling range 
/ is 2 {i.e., equal to the number of harmonics of the 
potential). 

Our purpose is to study the transport properties in 
this system. We limit ourselves to the case of a square- 
wave force switching alternatively between ±F. Then 
the rectified velocity is defined as, 



7(p)r = l(p)+F + 7<P>- 



(13) 



As we consider adiabatic conditions (u -> 0), the {p)±f 
are the corresponding stationary velocities; they can 
be calculated immediately with the general result (7) . 

In ref. [11] we considered the temperature depen- 
dence of the ratchet current. In this article, we study 
the force dependence of r y{p} 1 at fixed temperature. A 
high enough temperature is chosen to circumvent va- 
lidity problems of the master equation (see below) . 

4.1. Classical Limit 

First, it is convenient to understand the classical 
phenomenology. We plot in figure 1 the effect of finite 
inertia on the curves j{p)i vs. J 7 at a fixed temperature 
k B T/E = 1. 

The results obtained can be summarized as follows: 
(i) the rectified velocity is positive, i.e., the net drift is 
to the "easy" side; (ii) at high 7 the results converge 
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Fig. 1. Rectified velocity vs. force for a classical particle in 
a ratchet potential (a square wave ±F is used). Results are 
shown for various values of the damping 7, as well as the 
analytical overdamped result. The vertical lines mark the 
two critical forces in the deterministic case F^ [Eq. (14)]. 

to the known overdamped analytical result; (iii) finite 
inertia (7) shifts the curves to lower forces and narrows 
them; (iv) at low enough damping the net current drops 
to zero. 

There are two main points to be understood here: 
the dependence on 7 of the efficiency of the rectification 
and the window of forces where it occurs. Recall that 
in the zero noise limit (T = 0) , there exist two critical 
forces (for r > 1/4), 



\Fc/V \ 



1 + r, 



\F C + /V \ 



(14) 



where /Vo is the critical force for barrier dissapper- 
ance to the right and left sides. In our case r = 0.44, 
then F'/Va > F+/V . Thus in the T = limit we 
expect a current to the right (i.e. positive). 

Borromeo and co-workers studied the classical in- 
ertial, deterministic (T = 0) case [16]. They found 
that the net current starts at F^/Vb. Then rectifica- 
tion grows quickly with F, up to |F C - /Vb|. From this 
force the rectification starts to decrease with F, tend- 
ing to zero at large forces. They also observed that in- 
ertia narrows the curves, and fits them to the window 
(F c + /Vb, F~/Vq). In addition, for intermediate damp- 
ing the ratchet efficiency grows inversely proportional 
to the damping. Finally, at very low damping the rec- 
tified velocity drops to zero. Thus, our findings agree 
with the efficiency and narrowing characteristics of Ref . 
[16] for the noiseless limit. Nevertheless the tempera- 



ture smooths the curves in our case, allowing for recti- 
fied velocity below F+. 

This range below F c + is indeed quite rich already in 
the deterministic case. There are two more character- 
istic forces Fi(7) and F 2 (y) (Fl < F 2 < F c ). Below Fl 
the locked solution ((p) ~ 0) is the attractor globally 
stable. Above F2 the stable attractor is the running 
solution ((p) ~ F/7). Between F\ and F2 there exist 
bistability between these solutions. However, both for 
locked and running solutions there is little assymetry 
in the response [as then (p) is nearly independent of 
V(x)]. Therefore, the assymetry in the response con- 
centrates in the range Fi < F < F2. Then, since both 
Fi and F2 decrease with 7 [15, p. 330], the window with 
the maximuun rectification is shifted to lower F as 7 
is decreased, in accordance with our results. 

4.2. Quantum corrections 

After the exploration of the classical limit, we pro- 
ceed to make the system quantal, by decreasing K (see 
section 2), and see how this affects the ratchet effect. 
The physical meaning of K can be seen calculating the 
bands of the closed system, which depend only on K 
[11]. The number of bands grows with K (as a rule of 
thumb the number bands below the potential is ap- 
proximately K/2) while the band width increases with 
decreasing K. We plot in Fig. 2 the bands for the case 
K = 5. 
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Fig. 2. Energy bands for K = 5. The potential profile 
V[x) = — Vb [ sin x + (r/2) sin(2x)l is plotted (dashes) to 
show the number of bands below the barrier. In the calcu- 
lations we use r = 0.44 which smooths the potential profile. 
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We show results for the rectified current j(p) T in 
Fig. 3. We present curves at ksT/Eo = 1 and 7 = 
0.05 for different values of K. It is observed that the 
deviation from the classical case is systematic in K. 
First, we find a reduction of the ratchet effect at low 
forces. The reduction continues at intermediate forces, 
while at high forces an enhancement is observed, with 
respect to the classical case. Finally, a slight shift of 
the peaks of the curves to higher forces is found as K 
is decreased. 

The reduction/enhancement of the rectified velocity 
with respect to the classical limit was already investi- 
gated in [11]. These deviations can be understood in 
terms of quasiclassical corrections to classical trans- 
port, which is modified by tunnel events and over- 
barrier wave reflection (reflection of particles with 
energy above the barrier; those energies are attained 
by thermal activation) . Modelling with an asymmetric 
saw tooth potential one finds that tunnel events are 
more frequent in the hard direction that in the easy 
one for moderate to weak amplitude forces (leading to 
reduction of the ratchet effect)[l,ll]. On the contrary, 
the phenomenon of the wave reflection is less intense 
in the easy direction [11]. 

In particular, at F ~ 7 we obtained in [11] that, 
depending on the inertia and temperature, the reduc- 
tion/enhancement is a consequence of the "competi- 
tion" of tunnel events and wave reflection. Increasing 
the force both the contributions of wave reflection and 
tunnel transmission grow [11, App. F]. However wave 
reflection grows like ~ F, while tunnel transmission 
goes exponentially with F. Thus, when F > 7, but not 
too high, tunnel contributions are favored more than 
wave reflection. As a consequence, one finds a reduction 
of the rectified velocity in the range of inertia where 
an enhancement at F < 7 takes place. Finally, at high 
enough forces the barrier is sufficiently lowered. Then, 
there are many particles with energies above the bar- 
rier (experiencing wave reflection) . In addition, at high 
forces tunnel events also favor the enhancement [1] . 

As for the shift to higher forces, Shushin and Pollak 
showed [17] that quantum corrections produce a sup- 
pression of the diffusion in the underdamped regime. 
This effect yields an effective increased barrier, which 
manifest itself in a shift to higher forces in the quan- 
tum curves. 

In figure 3 (see the inset) we also observe a small cur- 
rent at F = 0. We believe that this "little" violation 



of the second law (transformation of thermal fluctua- 
tions in net current in absence of forcing) is a conse- 
quence of the approximations done in the derivation 
of the Caldeira-Leggett equation [8] (since the starting 
formalism is thermodynamically consistent). 
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Fig. 3. Rectified velocity vs. force for a quantum particle in 
a ratchet potential (excited with ±F). Results are shown for 
various values of the quantumncss parameter K (see main 
text). The damping and temperature are fixed at 7 = 0.05 
and ksT/Eo = 1. We show for reference the classical case 
(solid). Inset: Amplification of the very low force range. A 
small negative rectified current it is observed at F = for 
the quantum curves. 



5. PHASE SPACE REPRESENTATION 

We plot a (stationary) Wigner distribution in fig. 4. 
In order to see clearly quantum effects we choose the 
lowest K case consider (K oc Sq/K). 

The white islands in the graph represent zones with 
negative values of the Wigner function. Apart from this 
negative zones, the Wigner distribution lies mainly on 
its classical counterpart. 

Negative zones are a consequence of quantum in- 
terference. Thus, one would expect that decoherence 
yields positive Wigner functions. In fact, in the har- 
monic dissipative case, it was shown that after a finite 
time the Wigner function becomes positive [18] . How- 
ever, our plots show that more general potentials allow 
for non positive Wigner distributions at long times. 

By decoherence we mean that an intial reduced den- 
sity matrix, g(to) evolves with time to a diagonal one 
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all the structure of the stationary solutions obtained 
solving the Caldeira-Leggett equation for small 7/T. 




Fig. 4. Wigner function for a particle in a ratchet potential. 
The quantumncss parameter is K = 5, the damping is 
7 = 0.05 and F = 0. White "islands" correspond to zones 
of negative W (~ — 10~ 5 ). These negative zones arc of size 
comparable to those of Hamiltonian case [11]. Note the 
spatial asymmetry induced by the ratchet potential. 

(in a preferred basis) , when the system is coupled to the 
environment [19]. The basis {|n)} in which q becomes 
diagonal, Q = Y1 l c ™| 2 | n )( n l> it i s the so-called pointer 
basis [19]. The preferred basis depends on the Hamil- 
tonian of the system (Hs) and the interaction between 
the system and the bath. For example, in linear Quan- 
tum Brownian Motion (free particle and harmonic os- 
cillator) , the pointer states are localized in phase-space, 
the so-called coherent states (familiar from quantum 
optics). Mixtures of these states yield positivite W. 

The fact that our stationary W retains the negative 
islands (after decoherence) , implies that such coherent 
states cannot be the pointer basis in our case. Actually, 
in the limit of weak damping and high temperature, one 
expects that the stationary solution of (2) corresponds 
to the canonical density matrix, 

Q=^eT PHs = ^e- 0E *\k)(k\. (15) 

Accordingly, the pointer states would correspond to 
eigenstates of the system Hamiltonian. 

This can be additionally supported by the follow- 
ing. Consider the approximate Bloch functions that 
can be analytically obtained in the extreme quantum 
case (K < 1) [11]. Weighting the corresponding \k){k\ 
by the thermal factor e~ l3Ek and integrating over fc's 
in the lowest bands, one obtains Wigner functions with 



6. CONCLUSIONS 

The continued-fraction method is an established 
technique to solve classical Fokker-Planck equations. 
Admittedly, compared with direct simulations, the 
method has several shortcomings (it is quite specific 
of the concrete problem; the stability and convergence 
fails in some ranges of parameters, etc.). However, it 
also has valuable advantages (it is free from statistical 
errors, its nonperturbative character, high efficiency, 
etc..) [15]. 

When transferred to the quantum case the method 
inherits these shortcomings plus the problem of the 
critical election of the basis and recurrence index [11]. 
However, no quantum Langevin simulations are avail- 
able, while numerical solutions of master equations are 
computationally expensive. Notable exceptions are the 
different versions of the quantum stochastic calculus 
[20], which can be viewed as a sort of quantum coun- 
tenpart of the classical Langevin simulations. 

Let us emphasize that the continued-fraction tech- 
nique gives the exact solution of the quantum master 
equation ( Wigner distribution), allowing the calcula- 
tion of any observable [vd. Eq. (7)]. This method has 
also "spectral advantages" since it does not require 
the calculation of the eigenvalues, circumventing the 
demanding problem of continuum spectrum. Besides 
Quantum-Classical connection is attained in a natural 
way (tuned by K). Finally the method enjoys the ad- 
vantages of working in phase-space. 

In this work, we have studied the ratchet effect in 
the quantum case for high temperature/weak damping 
conditions. Starting from the classical case, we have 
proceed to make the system quantal. The deviations 
from the classical limit have been understood in terms 
of tunnel events and thermally activated over-barrier 
wave reflection. 

The representation of the Wigner functions has 
shown that, contrary the harmonic case, ratchet (co- 
sine) potentials do not exhibit positive Wigner func- 
tions at long times. We have argued that in the weak 
damping and high temperature regime the pointer 
basis are the eigenstates of the system Hamiltonian, 
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yielding in our case a non positive W. 

We finally mention that the limitations of the mas- 
ter equation considered do not allow us to study the 
low temperature regime, where quantum corrections 
will be more important. These problems are brought to 
the fore with an unphysical violation of the second law 
of Thermodynamics. This situation constitutes a mo- 
tivation to develop corrections to the Caldeira-Leggett 
equation, which could hopefully be solved with the con- 
tinued fraction method. 
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